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ABSTRACT 

We consider non-local effects which arise when radiation emitted at one radius of 
an accretion disk either heats or cools gas at other radii through Compton scattering. 
We discuss three situations: 

1. Radiation from the inner regions of an advection-dominated flow Compton 
cooling gas at intermediate radii and Compton heating gas at large radii. 

2. Soft radiation from an outer thin accretion disk Compton cooling a hot one- or 
two-temperature flow on the inside. 

3. Soft radiation from an inner thin accretion disk Compton cooling hot gas in a 
surrounding one-temperature flow. 

We describe how previous results are modified by these non-local interactions. 

We find that Compton heating or cooling of the gas by the radiation emitted in the 
inner regions of a hot flow is not important. Likewise, Compton cooling by the soft 
photons from an outer thin disk is negligible when the transition from a cold to a hot 
flow occurs at a radius greater than some minimum i?tr,min- However, if the hot flow 
terminates at i? < Rtr,mirn non-local cooling is so strong that the hot gas is cooled to a 
thin disk configuration in a runaway process. In the case of a thin disk surrounded by 
a hot one-temperature flow, we find that Compton cooling by soft radiation dominates 
over local cooling in the hot gas for M ^ and R ^ W'^Rschw As a 

result, the maximum accretion rate for which an advection-dominated one-temperature 
solution exists, decreases by a factor of ~ 10, compared to the value computed under 
an assumption of local energy balance. 


Subject headings: accretion, accretion disks - radiation mechanisms: nonthermal - 
radiative transfer 


1. Introduction 


Hot optically thin accretion flow solutions are often used to model X-ray binaries and active 
galactic nuclei. In addition to the original two-temperature model of Shapiro, Lightman & Eardley 
(1976) and the corona model of Haardt &: Maraschi (1991), two new classes of advection-dominated 
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models have been proposed recently: a two-temperature solution (Narayan & Yi 1995, Abramowicz 
et al. 1995), and a one-temperature solution (Esin et al. 1996). 

Advection-dominated models differ from the usual accretion solutions in that only a fraction 
of the viscously dissipated energy is radiated locally in the disk. The remainder of the generated 
energy is stored in the gas as entropy, resulting in a very hot, optically thin, quasispherical 
flow. If the accreting object is a black hole, as it is believed to be the case for several X-ray 
binaries and all active galactic nuclei, the stored energy is carried by the gas inside the horizon, 
so that the luminosity of the system can be orders of magnitude smaller than is predicted by 
standard accretion theory. The gas in advection-dominated solutions cools primarily through 
bremsstrahlung, synchrotron and inverse Compton scattering processes. Since the accreting 
electrons are heated to temperatures of up to 10^^ K, the spectra of such systems span ~ 10 
orders of magnitude in energy, from synchrotron radio emission to ~ 400/ceY X-rays produced via 
Compton cooling of hot electrons. Because two-temperature advection-dominated solutions have 
low radiative efficiency and hard power-law spectra, they have been successful in reproducing the 
observed properties of various low-luminosity X-ray sources including Sagittarius A* (Narayan, 
Yi, & Mahadevan 1995), the black hole nova A0620-00 in quiescence (Narayan, McClintock & Yi 

1996) , and the central source in NGC 4258 (Lasota et al. 1996). They also appear to be relevant 
to luminous systems (Narayan 1996). It is not yet clear if one-temperature solutions are relevant 
to any observed system. 

Though detailed numerical spectra of advection-dominated ffows include global radiative 
transfer effects (Narayan, McClintock & Yi 1996; Narayan 1996; Narayan, Barret & McClintock 

1997) , analytical calculations of the physical properties of these accretion solutions, which we 
briefly review in §2, generally ignore the effects of non-local radiative transfer (Narayan & Yi 1995; 
Esin et al. 1996). This is a dangerous approximation. The optical depth for electron scattering in 
such flows is generally much lower than unity, and therefore, radiation emitted at one radius can 
potentially change the energy balance of the gas at quite different radii. It is the purpose of this 
paper to investigate the importance of this effect. We analyze three different but related issues. 

Eirst, in the advection-dominated models the bulk of the luminosity is emitted from the 
innermost region of the disk {R < W0Rschw)i where most of the gravitational energy is dissipated 
and the cooling efficiency of the gas is higher. The radiative flux from this region can have two 
effects on the rest of the flow, as we discuss in §3.1. The hot radiation can Compton heat (§3.1.1) 
the outer regions {R > WOORschw), where the gas is cooler (Shvartsman 1971, Ostriker at al. 
1976, Grindlay 1978). On the other hand, cold but numerous synchrotron photons can contribute 
to cooling of the gas (§3.1.2) closer in {lOORschw < R < lOOORschw)- 

Next, the hot solutions based on the one-temperature advection model (Esin et al. 1996) do 
not extend down to the last stable orbit for M > At accretion rates higher than 

this limit the only possible configuration is a standard thin disk (Shakura & Sunyaev 1973; Frank, 
King & Raine 1992), surrounded by a hot flow. Since most of the energy is dissipated in the thin 
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disk, the number density of the cold photons emitted in the inner regions can be much greater 
than the number density of photons emitted locally by the hot gas, and we expect that the hot 
flow will suffer significant extra cooling via Compton scattering of the thin disk photons. This 
issue is discussed in §3.2. 

Finally, both one- and two-temperature advection-dominated solutions can be surrounded by 
a thin accretion disk, with the boundary between them at some radius Rtr, where the value of Rtr 
depends on the parameters of the model. Here again the cold photons emitted by the thin disk 
can Compton cool electrons in the nearby hot flow. We analyze this situation in §3.3. 

Just as the radiation from the thin disk affects the energy balance of the hot flow, the cold gas 
itself is affected by irradiation from the advection-dominated zone. In §3.4 we investigate whether 
this process affects our conclusions from §3.2 and §3.3. We summarize and discuss the results in 
§4. 


2. Advection-Dominated Accretion Flows 

In this section we review the basic equations, derived by Narayan & Yi (1995), that describe 
local properties of advection-dominated accretion flows in terms of five basic parameters: mass 
of the accreting object, M, mass accretion rate, M, radial distance from the accreting object, R, 
standard viscosity coefficient, a (Shakura & Sunyaev 1973), and the ratio of gas pressure to total 
pressure, f3. 


2.1. Self-Similar Solution 


The total pressure in a vertically averaged, axisymmetric accretion flow is given by 

P = pel, 


( 2 - 1 ) 


where p is the gas density and Cs is the local isothermal sound speed. We assume that the gas 
is in equipartition with a tangled magnetic field, so the total pressure is the sum of the gas and 
magnetic pressures, p = Pg + Pm- In the ionized gas, the expression for the gas pressure is 


Pg = Ppcl = 


pkTi pkR 


+ 


( 2 - 2 ) 


PilTlu Perriu 

where Tj, Tg, and pi, pe are temperatures and molecular weights of ions and electrons respectively. 
The magnetic pressure due to the isotropic tangled magnetic field is given by 


= (1 - (5)pc, = 


247r 


(2-3) 
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Note that this equation differs from the one adopted by Narayan &: Yi (1995) by a factor of 1/3. 

Narayan & Yi (1995) have derived a self-similar solution for advection-dominated accretion 
flows that obey mass, energy, and momentum conservation. Combining this solution with the 
equation for the gas pressure, and scaling mass in units of the solar mass, M = tti-Mq, accretion 
rate in Eddington units, M = mMEdd = ^ 1-39 x lO^^m gs“^, and radii in Schwarzschild units, 
R = rRschw = T 2.95 x lO^m cm, we obtain 


1Q TTl q 

Ue = 1.33 X 10^^ -;-Cm 

(2-4) 

qt = 1.84 X 10^^ - erg cm ^s , 

(2-5) 

Ti + l.OSTe = 6.66 X 10^2— K, 

(2-6) 


where Ue is the electron density, is the rate of viscous energy dissipation per unit volume, and 
we have taken //j = 1.23 and = 1.14, as required for a gas composition of 75% H and 25% He. 
The parameters C 3 and e' are defined as 


C 3 = 


2(5 + 2e') 


1 + 


18 a^ 


(5 + 2e')2 


1/2 


- 1 


(2-7) 


and 

, _ 1 

' /17-1 y / ’ 


( 2 - 8 ) 


where / is the fraction of viscously dissipated energy that is stored in the gas, and the expression 
for the adiabatic exponent 7 is derived in Appendix A. 


2.2. Heating and Cooling Processes and Local Energy Balance 


2.2.1. Heating of Electrons 


Since ions are considerably more massive than electrons, viscous heating affects primarily the 
ions (Shapiro, Lightman, Eardley 1976; Phinney 1981; Rees et al. 1982). However, cooling of 
the hot plasma occurs mainly through electrons, so the energy has to be transferred in some way 
from ions to electrons. In their analysis, Narayan & Yi (1995) (following Shapiro, Lightman, &: 
Eardley [1976]) assumed that the energy transfer can occur only through Coulomb interactions at 
a volume transfer rate, qie, the expression for which was derived by Stepney & Guilbert 1983. Since 
qf (X rh (X Hi (Eq. |2-5| ) and qie oc UiUe (Stepney &: Guilbert 1983), at low densities. Coulomb 
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interactions are less efficient than viscous energy dissipation, giving rise to two-temperature 
accretion flows with ions much hotter than electrons. 

In addition to Coulomb scatterings, some non-thermal coupling mechanisms have been 
proposed recently (e.g. Begelman & Chiueh 1988). If such mechanisms are more efficient than 
Coulomb coupling, they can restore the thermal equilibrium between the ions and electrons, 
creating a one-temperature accretion flow (Esin et al. 1996). In that case, the question of energy 
transfer from ions to electrons is settled by simply setting Tj = Tg at all times. 


2.2.2. Local Cooling Processes 

The fraction (1 — /) of the energy generated by viscous heating must be radiated by the gas. 
Cooling of the hot plasma occurs via three processes: bremsstrahlung, synchrotron radiation, and 
inverse Compton scattering of low energy bremsstrahlung and synchrotron photons; 

Q Qbr ^sync T 

Bremsstrahlung and Line Emission The free-free emission in a plasma is produced in 
electron-electron and electron-ion interactions. The corresponding cooling rates were derived by 
Stepney & Guilbert (1983) and Svensson (1982). The total bremsstrahlung cooling rate per unit 
volume is the sum of the two contributions, q^r ~ ^ee Qei- 

At large radii, the gas is cold enough to allow line cooling, which becomes more efficient than 
free-free emission below 10® K. At these temperatures electron-electron bremsstrahlung emission 
is negligible, and we use a cooling function that includes electron-ion bremsstrahlung and line 
cooling (updated version; original version was described by Raymond, Cox, & Smith [1976]). 


Synchrotron Emission In a hot plasma with equipartition magnetic field, synchrotron 
radiation is an important cooling mechanism. An expression for the spectrum of optically thin 
synchrotron emission from a relativistic Maxwellian distribution of electrons was derived by 
Mahadevan, Narayan, & Yi (1996). However, below some critical frequency, z^c, the emission 
becomes self-absorbed and the optically thin emission formalism no longer applies. To estimate Vc 
we equate the synchrotron emission from a sphere of radius R to the Rayleigh-Jeans blackbody 
emission from the surface of that sphere. Below Uc the emission is optically thick and we use the 
blackbody estimate, while above Vc we use the expression derived by Mahadevan et al. (1996). 
Integration over frequency yields the total cooling rate per unit volume (Esin et al. 1996). 


Inverse Compton Scattering Both bremsstrahlung and synchrotron emission mechanisms 
produce mainly soft photons, which are then upscattered by the hot electrons in the accreting gas. 
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The optical depth to electron scattering is simply Te = 2ne(TTH, where H is the scale height of the 
flow determined through the conditions for vertical hydrostatic equilibrium. After each scattering, 
the energy of a photon is increased by a factor A = 1 + 46^ + 166^ (Rybicki & Lightman 1979), 
where 9e = /cTe/rUeC^ is the dimensionless electron temperature. The average ratio of the photon 
energy at escape to its initial energy is given by (Esin et al. 1996) 

f] = ^^^[1 — P{jm 4-4)] + rjmaxPUm 4-1, s), (2-10) 

where P{a,x) = [l/r(a)] Jq is the incomplete gamma function and 

s = Te + rl, (2-11) 

jm = ln(r/max)/ln(A), (2-12) 

Vmax = SkTe/hu. (2-13) 

The synchrotron spectrum is strongly peaked at Uc, so we can estimate the Compton cooling 
rate due to the scattering of synchrotron photons as q'^ ^ = (r?(^'c) — l)<?^nc- On the other hand, the 
spectrum of bremsstrahlung radiation is practically flat between the synchrotron self-absorption 
frequency I'c (we found that the bremsstrahlung self-absorption frequency is generally lower than 
Uc for the range of accretion flow parameters that is relevant for this work) and exponential cut-off 
at kT/h. Below Uc free-free emission falls off as Rayleigh-Jeans blackbody and Comptonization is 
not important, since the photons are more likely to be absorbed than scattered. Consequently, we 
can estimate the emission per unit frequency ss ^ cs . Then the Compton cooling 

- kT 

rate is given by the integral q^^^^\q{v) - l]dv. 

The total Compton cooling rate is simply q^ = q^s 4- <lcbr- 


2.2.3. Energy Balance 


The local energy balance in the gas requires that a fraction (1 — /) of the viscously dissipated 
energy is radiated away. This gives us the condition 

= (2-14) 

which together with the equation of state and self-similar solution described in §2.1 is sufficient 
to obtain a closed set, necessary to solve for all one-temperature flow parameters. For a 
two-temperature solution, we need another equation to determine Tj and Eg separately. We obtain 
it by requiring that the net cooling rate has to be equal to the heating rate of the electrons; 


q = qie- 


( 2 - 15 ) 



2.3. Physical Properties of One- and Two-Temperature Flows 


The set of equations described in §2.1-2.2 allows all the major parameters describing the 
accreting gas to be solved for self-consistently for given values of M, M, R, a, and /3. The 
resulting flows are quasispherical, with H ^ R, optically thin, and quite hot. Their properties 
were described in detail by Narayan & Yi (1995) and Esin et al. (1996), and here we give a brief 
summary. 

As was shown first by Rees et al. (1982) advection-dominated accretion solutions are limited 
to values of rh that lie below a maximum mass accretion rate, mcrit{f')- Above this limiting value, 
the radiative efficiency of the gas is so high due to the increased plasma density, that the flow cools 
down to the standard thin accretion disk configuration. Figure 1(a) shows a comparison of the 
critical accretion rate as a function of R computed for one- and two-temperature models. At large 
radii, r > 10^, the two calculations give identical results, with rhcrit oc Note that previous 

papers by Narayan & Yi (1995) and Abramowicz et al. (1995) have merit oc in the outer 
regions. Our much stronger dependance of the critical accretion rate on radius is due to including 
line cooling, which is significantly more efficient than bremsstrahlung at low temperatures. Closer 
to the accreting object, merit decreases steeply as rherit oc in a one-temperature model, but 
is roughly independent of r in a two-temperature model (the slight decrease in merit seen on the 
figure is most likely an artifact of the self-similar assumption we used to derive the solution; when 
proper boundary conditions are used at the inner edge of the disk, this decrease is no longer 
present [c.f. Narayan 1996]). 

To understand why the two solutions behave in this way, we show in Figure 1(b) the radial 
temperature profile of the accreting gas for a IOMq black hole, with M = a = 0.1 

and (3 = 0.5. At large radii, the ion and electron temperature in a two-temperature flow are 
equal to each other and to the gas temperature in a one-temperature flow. In that region, the 
Coulomb energy transfer from ions to electrons is efficient and the two solutions are identical. 
However, in the region inside lOOORs'c/ito) the electron temperature in a two-temperature solution 
remains roughly constant, whereas the ions stay at near virial temperatures with Tj oc r~^. The 
gas in a one-temperature flow, on the other hand, necessarily finds an equilibrium temperature 
with Te = Ti = Tq. The temperature Tq obviously lies between Tg and R of the two-temperature 
solution. At ~ 10^*^ K, the rate of synchrotron and Compton cooling of the gas is very sensitive to 
the temperature of the electrons, and since electrons in one-temperature solutions are hotter than 
electrons in two-temperature solutions, the former has a much higher radiative efficiency. This in 
turn leads to lower merit- 

The difference in temperatures between the two solutions also leads to different values of the 
advection parameter /. Figures 1(c) and 1(d) show contours of constant / as a function of radius 
and accretion rate, for two- and one-temperature models respectively. Two-temperature flows 
are advection-dominated, / ~ 1 everywhere, except for a narrow region near the Merit boundary 
(shown as a heavy line). On the other hand, the inner regions of single-temperature accretion disks 



are cooling-dominated, / ^ 0.5, for rh ^ because they have higher electron temperatures 

In fact, one-temperature hot accretion flows at high M have radiative efficiency close to 15%, 
which means they radiate away nearly all the gravitational energy lost in the accretion process 
(Esin et al. 1995). By comparison, the efficiency of two-temperature accretion solutions can be 
orders of magnitude less. 


3. Effects of Nonlocal Radiation 


The formalism described in the previous section ignores the effects of non-local radiative 
transfer. It assumes that the photons emitted at one radius interact only with the ‘local’ gas, i.e. 
gas within a region of unit radial width in log (r). We can check the validity of this assumption by 
estimating the optical depth for electron scattering in the radial direction predicted by advection 
models. The optical depth from radius to r is simply Tr{rmin,r) = UeCTTRschwdr. After 
substituting the expression for Ue given by Eq. ( p-4| ) , where we take C3 ~ 0.3 (see Appendix B), 
and integrating, we obtain 


Tr{rn 


Trij’min) — 1'8 





(3-1) 


assuming that r » Tmin- Thus, two-temperature accretion disks are optically thin if 
M ^ 3 X lO^^aMgdd, which means that the photons can travel freely through the disk and be 
scattered by electrons far from where they were emitted, heating or cooling the gas as a result of 
scattering. 


One-temperature accretion disks are geometrically thinner in the inner region than two- 
temperature disks, because they are able to cool more efficiently (Esin et al. 1996). As a result, 
both He and are higher. However, even one-temperature flows are optically transparent below 
M ~ 3 X 


In the remainder of this section we attempt to quantify how Compton scattering of non-local 
radiation will affect the flow parameters calculated based on local analysis. In each case we 
estimate analytically the region in the parameter space where non-local effects can be important, 
then we compute numerically how the properties of the solutions are affected in this region. 


3.1. Photons from the Hot Inner Region 


Both one- and two-temperature accretion flows produce most of their emission from within 
lOOR^cW) where most of the gravitational potential energy is released. The inner region also 
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has the highest electron temperatures so that the emerging photons are hotter than the gas in 
the surrounding flow. When this radiation interacts with the gas in the outer regions of the disk, 
two processes can take place. Where the gas is relatively hot, less energetic photons will cool the 
electrons through inverse Compton scattering. On the other hand, where the temperature of the 
electrons is less than the effective temperature of the radiation, the electrons will be heated via 
Compton scattering. 

The amount of the heating or cooling can be easily estimated if we assume that each photon 
undergoes at most one scattering at any given radius (a reasonable assumption, since the Compton 
y-parameter (Rybicki & Lightman 1979) of the accreting gas at i? > lOORsc/iio is always less than 
unity). With this approximation, the rate of energy transfer from photons to electrons per unit 
volume, via Compton scattering, is just 

Qni = crTcnee~'^^ j AE,y dv ergs“^cm“^, (3-2) 

where is the number density of the incident photons per unit frequency and AEy = 
is the energy lost by the photon in each scattering. When q^i is positive, the electrons are heated 
by the incident photon flux, when q^i is negative, the electrons are cooled. We approximate n^^y 
as riry^y = Ey jcku ~ Ly /E}ckv, where Ly is the luminosity per unit frequency emerging from 
the region within 100i?5c/i^. AEy is given by Rybicki & Lightman (1979) 

AEy = E^ iri - -E'7,fin - - AkTe). (3-3) 

rrieC^ 

Since AEy can be at most as large as the original photon energy, Eq. ( ^-3| ) is valid only for 
nonrelativistic photons with hv < mgC^- To extend this formula to higher photon energies, we 
have to replace hvlnif^c? by hv/{me(? — hv), to ensure that AEy < hv. However, numerical 
calculations of qni (described in §3.1.1 and §3.1.2) show that this modification has a very small 
effect on the final result and so for simplicity we use Eq. (p-3|) in our analytic estimates. 


3.1.1. Heating of the Gas 


Let us consider first the rate of heating of the gas by the energetic photons. Erom Eq. ( p-3| ) 
it is clear that only photons with energies hn > 4:kTf>{r) can contribute to the heating, whereas 
photons with lower energies will cool the gas. The high energy cutoff of the spectrum is determined 
by the electron temperature in the inner regions, Te^in, since inverse Compton scattering cannot 
produce photons with energies higher than ~ ‘ikTf.^in (Rybicki & Lightman 1979). Therefore, to 
obtain an upper limit on the amount of heating, the integral in Eq. (^) has to be evaluated from 
v = AkTf,{r)/h to v = 3kT(>^in/h. Clearly then, significant amount of Compton heating can occur 
only if Te(r) <C Te.in, which restricts us to the region outside lOOORschw (see Figure l[b]). 
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The precise shape of depends sensitively on the value of the accretion rate. However, a 
reasonable upper limit to the effect of heating can be obtained by assuming that the spectrum is 
flat, i.e. equal energy is emitted per logarithmic interval in frequency, and that the luminosity 
emitted from the inner region is approximately equal to the total luminosity of the disk. Since the 
spectra generally span ~ 10 orders of magnitude in frequency (e.g. Narayan, Yi, &: Mahadevan 
1995; Lasota et al. 1996), we can write Ly ~ L/(101oggl0i^), with L being the total emitted 
luminosity. Using this approximation, and substituting Eqs. (|2-4D and (|3-3|) into Eq. (^-2| ) , we 
obtain an approximate expression for the rate of heating caused by non-local radiation from 
smaller radii: 


where we assumed C3 
result. 


qti 1.1 X 10 


-18 


me 




m^ar'^l'^ V10® K 


- e,2n 


— 1 — 
erg s cm 


(3-4) 


~ 0.3 (Appendix B) and used the fact that Tg to simplify the final 


Compton heating of electrons by hot photons becomes important when it is comparable to the 
rate of energy transfer from ions to electrons, so to find the region of the parameter space where 
non-local heating must be taken into account, we impose the condition q'^i/qte = Qni/Q~ — 1- The 
expression for q~ is derived in Appendix B, and the left-hand side is given by Eq. (p-4|) . 


Eor two-temperature accretion disks, this condition requires that 

ioooy/2 /io9k\ ^ 
r J ” 


m ^ 2.3 X 10 


-2 


a 

or 


e,2n 


(3-5) 


for r > 1000. The values of m and r for which this relation is satisfied lie above the dashed line on 
Eigure 2(a), computed assuming a = 0.3, (3 = 0.5 and Tr <C 1. One can see that the region where 
significant Compton heating is expected to occur is small. 


In reality, since we used upper limits in estimating our analysis overestimates the 
magnitude of this effect. In order to find the exact values, we have solved the complete set 
of equations discussed in §2 numerically, and computed the values of electron density, plasma 
temperatures, and magnetic field strength in a two-temperature accretion flow for various values 
of r and m. Using these parameters we calculate the local cooling rate, q~, and the shape of the 
spectrum of radiation emerging from the inner region (see Narayan, Barret &: McClintock [1997] 
for details). Then the integral in Eq. (3^2) can be computed numerically using exact values for 
Tie, and AE, to obtain q'^^. On Eigure 2(a) we have plotted the contours of constant 

q^i/q~ as dotted lines. It is clear from this figure that non-local radiation mainly Compton cools 
the accreting gas, since the values of this ratio are for the most part negative. Where Compton 
heating does occur, it makes up at most a few percent of the energy transferred to electrons from 
ions, and can therefore be completely neglected. 


Eor one-temperature flows, the total disk luminosity can be estimated from Esin et al. (1996) 

1.3 X 10^® mm erg s“^ for m^l0“^a^. 


Li-t ^ 


1 


3 ^ ^q 43 g^gg -1 


for m ^ 10 ^o?. 


(3-6) 
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Then the condition q^i > q requires that 


„ a 4.6 MO- 



(3-7) 


for r > 1000. The region where Eq. (3-7) is valid lies above the dashed line on Figure 2(b), 
for a = 0.1, (3 = 0.5, and <C 1. Note that the area affected by non-local heating is larger 
than the similar estimate for a two-temperature flow on Figure 2(a). This happens because 
one-temperature disks are more luminous than their two-temperature counterparts. However, this 
result is again an overestimate. As we did in the case of a two-temperature flow, we computed 
the ratio q -\-ni /<?” numerically, solving the flow equations appropriate for a one-temperature 
plasma. The constant value contours of this ratio are shown on Figure 2(b) (dotted lines). The 
results show that Compton heating is indeed stronger here than it was in a two-temperature gas. 
It is most important for values of the mass accretion rate m ~ nicriti'r = 3), since at higher m 
the hot flow can not extend all the way to SRschw: and at lower rh the optical depth is low and 
Comptonization is inefficient, so that fewer high energy photons are produced, capable of heating 
the gas in the outer flow. However even at its maximum, non-local Compton heating does not 
exceed ~ 10% of the local heating of electrons via viscous dissipation. 

We conclude that Compton heating can be neglected within the advection-dominated zone 
both for two- and one-temperature flows. 


3.1.2. Cooling of the Gas 


The photons with energies hn < AkT^, cool the gas through inverse Compton scattering. 
An upper limit on the amount of cooling induced by non local radiation can be obtained 


by assuming that all photons are emitted at <C AkTg/h, so that Eq. (3-3) reduces to 
AE ~ —hn{A9e + 160g) = hu{l — A), where the term 160g is added to account for relativistic 
effects at 0g ^ 1. Since at r > 100 the electron temperature in both one- and two-temperature 
flows varies as Tg ~ 10^^/3K/r (see Figure 1(a) and Eq. |^]), the quantity A — 1 can be written 
as A — 1 ~ 670P/r A (670/3/r)^. With this simplification, the integral in Eq. (^) is trivial and 
the expression for the cooling rate evaluates to 


Qni = 5-0 X 10 


-17 


670/3 / 670/3 \ 2' 

r \ r J 


me 


"L 


m^ar 


7/2 


erg s ^ cm ^, 


(3-8) 


where we have set = \qni\- 


Substituting Eq. (|B5|) for L and imposing the requirement q^Jq > 1, we find that in 
two-temperature flows photons from the inner region will contribute to the cooling of the gas at 
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r > 100 when 


m ^ 6 X 10 ^ 





1 + 3.4 




(3-9) 


For a = 0.3, /3 = 0.5, and <C 1, this inequality defines a region that lies above the dashed 
curve on Figure 2(c). We see that the extra cooling can be important only for high values of 
the accretion rate, when the electron density in the flow is large enough to scatter the outgoing 
photons from smaller radii. 


The results of an exact numerical calculation are also shown on Figure 2(c). We have 
computed q~i and q~ as described in §3.1.1, and plotted the curve with q~i/q~ = 1 as a solid 
line. The shaded region inside this curve, where non-local cooling can affect the accretion flow, 
is very small. To investigate how the flow properties change in the affected region, we have 
included q~i into our numerical calculations as an extra cooling mechanism, replacing Eq. (2-9) by 
(1~ = %r Q^ync + + Qni' ^ result, the equilibrium electron temperature of the accreting 

gas and the advection parameter / decrease slightly, to compensate for the increased cooling. 
However, since most of the emission produced by the flow comes from within lOOi^scW; there is 
practically no change in the observed spectrum. 


To estimate where non-local cooling is significant in one-temperature flows, we substitute Eq. 
(3-6) into Eq. (|3-8D and compare the result to the local cooling rate, q~, obtaining the following 
condition on m and r 


m ^ 1.2 X 10 




-1 


(3-10) 


for r > 100. The region of the parameter space where this inequality is satisfied lies above the 
dashed line on Figure 2(d), computed for a = 0.1, j3 = 0.5, and 1. The results of a detailed 
numerical computation are shown as a shaded region on the same figure. 

From comparison of Figures 2(c) and 2(d), it is clear that non-local cooling is much more 
important in one-temperature than in two-temperature accretion flows. However, in both cases, 
the only effect is to decrease the gas temperature and the advection parameter in the affected 
regions of the flow, by a modest amount. However, these regions do not produce most of the 
observed emission. Therefore, we conclude that the radiation emitted in the inner regions of a hot 
flow does not have significant observable effect on the gas in the outer regions. 


3.2. Cold Photons from the Inner Thin Disk 


In a one-temperature accretion flow model, the value of the critical accretion rate, merit■, 
increases as oc for r 1000 (see Fig. l[a]). For m > nicrit, the only equilibrium solution 
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is the thin accretion disk. As a result, when the value of the accretion rate is in the range 
^ rh ^ the only possible flow configuration is a thin disk surrounded by a hot 

quasi-spherical flow. Cold photons emitted by the thin disk will Compton cool the hot gas around 
it. This is a similar problem to the one we discussed in §3.1.2. However, since the luminosity 
considered here is higher because of the higher efficiency of the thin disk, we expect that non-local 
Compton cooling will be more important. 


In this scenario, non-local cooling is no longer limited to the outer regions of the disk, where 
the accreting gas has low density and temperature and therefore, multiple scatterings of the 
incident photons are not important. In fact, when we compute the Compton y-parameter we 
find that it can be much greater than unity for r <C 100. Thus, to estimate the rate of non-local 
Compton cooling, we need to modifying the expression for AE^ to include the effect of multiple 
scatterings: 

AE^ = E^^in - ^7,fin = hu{l - r]{u)A), (3-11) 


with rj given by Eq. ( ^-10 ). In general, the value of rj depends on the frequency of the scattered 
radiation, but for simplicity we set it equal to some average value fj. As in §3.1.2, we assume that 
the radiation is cold and the factor hi' is small as compared to kT^. With that approximation, the 
Compton cooling rate becomes 


^nl = crmeC '^^Fc{riA - 1 ). 


(3-12) 


where Fc = is the cold photon flux and Lc is the total luminosity of the thin disk. Since 

the standard thin disk has approximately Keplerian energetics, we can estimate Lc as the energy 
loss between the outer and inner radii of the disk: 


GMM 

‘^kischw 




(3-13) 


Combining Eqs. (|3-12D and (|3-13|) , we obtain the following expression for the cooling rate 


Qnl ^ 1-5 X 1023 




(^^-1) Q 



(3-14) 


Compton cooling becomes important when it exceeds local cooling of the gas. Comparing Eqs. 
( 3-14 ) and (|B^ ) we find that is greater than q at all radii, provided that rh ^ 10 
(necessary for the thin disk to be present) and fjA — 1 > 1. Thus, if the thin disk is present, it will 
always Compton cool the hot flow around it. 


To obtain more detailed results, we computed the quantities q~ and q~i numerically. In our 
calculation, we assumed that all the radiation from the thin disk is emitted at the frequency 
corresponding to a blackbody of temperature Tbb — {3GMM/SirR^aB)^^^ (Frank, King, &: Raine 
1992) evaluated at R = 3Rschw, and computed the corresponding fj at each radius. The region 
where q^i/q~ > 1 is shaded on Figure 3(a). The flow in this region will be cooled significantly 
by photons from the thin disk. Including non-local Compton cooling directly in our calculations 
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shows that because of the increase in the cooling rate, the gas near the inner edge of the hot flow 
collapses to a thin disk, and the critical accretion rate line itself moves downwards. This means 
that the thin disk extends slightly further outward and produces more cold radiation. After a 
few iterations, this process converges; the new Merit is shown on Figure 3(a) as a dashed line. 
Due to extra cooling, the maximum accretion rate for which hot one-temperature solution exists, 
decreased by a factor of ~ 10. For M < Merits where hot solutions are still allowed, the equilibrium 
temperature of the accreting gas decreases (Figure 3[b]), so that the observed spectrum contains 
less high energy flux. 


3.3. Cold Photons from the Outer Thin Disk 


Outside 1000ii5c/i«;, the critical accretion rate for both one- and two-temperature solutions 
decreases with radius as merit oc . Therefore, for any value of m there is a critical radius 
'^'crit beyond which the only possible flow configuration is the thin disk. Therefore, we expect 
to have accretion flows where the gas forms a cool thin disk for r > Teriti but switches to a hot 
advection-dominated state for r < Terit- Many of the models published in the literature (e.g. 
Narayan et al. 1996) are of this form. 

In such configurations, the photons produced in the cold outer flow may be expected to 
Compton cool the hot gas inside Verit- The cooling rate per unit volume is given by Eq. ( p-12 ), 
the only unknown being the cold photon flux, Fe- 

In general, Fe is a function of both radius, i?, and the height above the plane of the thin disk, 
z = RsinO. Taking into account the full geometry of the problem, we write Fe{R,z) as a double 
integral: 

Fe{R, z) = r dR' 2 r R' d(j) 
jRtr Jo 

where Rtr is the transition radius between the hot and cold parts of the accretion disk, Tc{R') is 
the surface temperature of the thin disk at radial position R', x is the distance from the emitting 
element in the thin disk to the scattering element in the hot flow, cj) is the azimuthal angle, and 
sin^ = zjx. Since we are working with a height-averaged set of equations, we replace Fe{R,z) 
with its average over a spherical shell of radius R: 

Fe{R)= cosOdO dR'2 R'd(t> ^ ^ (3-16) 

Jo JRtr Jo TTX"' 


The emission per unit area of the thin disk can be approximated as (TbTc — ^GMM/SttR^ 
(Frank, King, & Raine 1992), and substituting expressions for x and '0 in terms of R, R' and 9 we 
obtain 


Fe[R) = 


3GMMR 


' Rti 


dR' 


fK/2 


d6 cos 9 sin 9 


dej) 


(R'2 + _ 2 R'R cos 9 cos 0)3/2 


(3-17) 


0 
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This integral has to be evaluated numerically, but a reasonable estimate can be obtained by 
replacing sin if; by sin'll, where iIj is some effective angle averaged over the entire spherical shell. 
Since sin ijj = RsinO/{R'^ + — 2i?'iicos 6 cos setting (j) = 0, 6 = cos“^ R/Rtr, and R' = Rtr 

gives an upper limit of sin ip = R/Rtr- Then the three integrals can be evaluated analytically: 


Fc{R) 


3GMM R 
dvr Rtr 


1 


2RtrB? 


+ 


d 2 p2 
IL '^tV 

m/rR^ 


In 


ritr 


— In 



(3-18) 


In the limit when R ~ Rtr, Eq. ( p-lSj ) reduces to Fc ~ 3GMM/SttR/j-R] when R <C Rtr, 
Fc ~ GMMR/ATrRfr, so that when R/Rtr 0, the cold photon flux becomes completely 
negligible. 


Clearly, Compton cooling will be strongest near the boundary between the hot and cold flows, 
so an upper limit on the amount of extra cooling can be obtained by replacing Eq. ( 3-18| ) by the 
limiting expression for Fc when R ~ Rtr- Substituting this result into Eq. ( 3-12| ) and imposing 
the familiar constraint q~i/q~ > 1, we find that external cooling is significant in the hot flow near 
the boundary if 

^ r y / a /0.5\2 - 1' " 


rtr ^ 200 


rtrJ VO.l 


P ) \ 


(3-19) 


for m ^ 10' 


and 


rtr ^ 10^ 




0.5 



{fjA{r) - if, 


(3-20) 


for m ^ 10 Since the critical outer radius for a hot flow is always greater than 10^ Rschw (see 
Eigure l[b]), this result shows that the thin disk radiation has no effect on the hot gas, so long as 
Rtr ~ Rcrit' 


The cooling from the outer disk can, however, be significant if Rtr Rcrit- There is no 
well-developed theory at this time for the transition radius Rtr and one could in principle imagine 
a situation where the gas remains in a cold state down to a much smaller radius than Rcrit- If we 
consider such a model, the extra cooling of the hot gas at i? ~ Rtr can be quite important. To 
show this, we computed the ratio q~i/q~ numerically, using the full expression for Fc (as given by 
Eq. | 3-18|| ). The results for one- and two-temperature hot flows are shown on Figures 4(a) and 
4(b) respectively, where we have plotted q~i/q~ evaluated at i? = Rtr for different values of m 
and rtr- We can see that external cooling begins to dominate when rtr 10^ - 10^, in reasonable 
agreement with our analytic approximation. 


As the transition radius decreases further, the ratio of non-local to local cooling at the 
boundary between a two-temperature hot flow and a thin disk becomes more extreme, increasing 
up to ~ 200 for high accretion rates, and causing a significant decrease in the equilibrium 
temperature of the hot gas. Note however, that only a relatively thin layer near the boundary is 
affected by extra cooling. Away from the boundary, the cold photon flux decreases linearly with r 
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and the ratio of non-local to local cooling falls steeply with decreasing radius (see Eqs. [p-12| and 

(3-21) 


Q 


,-4 


Because of this strong dependence on r, the net effect on the hot flow is negligible, as long as 
the hot gas in the boundary layer is in stable thermal equilibrium. However, the hot equilibrium 
configuration cannot be maintained for an arbitrarily strong cold photon flux. At some point, 
the total cooling in the boundary layer becomes so strong, that the gas temperature is forced 
significantly below the virial value, i.e. the hot quasispherical flow collapses into a cold thin disk. 
When this occurs, the transition radius effectively moves inwards, which means that the hot gas 
in the new boundary layer will experience even stronger external cooling, and will also inevitably 
collapse. In this fashion, a runaway instability develops that causes the entire hot flow to cool 
down to the standard Shakura & Sunyaev disk. 


We have computed numerically the values of the transition radius at which this instability 
occurs for different mass accretion rates. This was done by incorporating q~i directly into the 
calculations and decreasing rtr until the hot stable accretion flow solution at r ~ rtr disappeared. 
The resulting critical values of rtr are shown on Figure 5 as dashed lines. We see that this 
instability is an important consideration only for large values of the mass accretion rate. For 
M ^ 0.1 X a^MEdd, the thin disk has virtually no effect on the hot flow, irrespectively of where 
the transition between the cold and hot phases occurs. 


In a single-temperature hot flow, the value of the ratio qni/q~ evaluated at r = rtr never 
becomes much higher than unity (see Figure 4[b]), even for small values of the transition radius. 
The reason for this is that at r < 100, one-temperature flows are cooling-dominated, i.e. local 
cooling is very efficient. At r < rtr the cooling ratio behaves roughly in accordance with Eq. 
(3-21), so that external cooling is completely unimportant there. 


3.4. Irradiation of the Thin Disk by the Hot Flow 


In calculating the energy flux from the thin disk in §3.2 and §3.3 we used the standard 
(Frank, King, & Raine 1992) formalism and ignored the effects of irradiation of the cold gas by 
the advection flow. In reality, most of the hot photons incident on the thin disk are reprocessed 
there, so that the net emission from the thin disk is a blackbody at a higher temperature than 
predicted by the standard theory. To test whether irradiation is important, we need to compute 
the ratio of the incident energy flux to the viscous dissipation rate per unit area of the thin disk. 

Using cylindrical geometry for the hot flow, the incident energy per unit area of an inner thin 







- 17 - 


disk can be written as a triple integral 

poo 

Qinc{R') = dR 2 

J Rtr ' 


Jo 4:TTX^ 


(3-22) 


where R, R', z, x, and (j) are defined as in §3.3. Since the thin disk is surrounded with a 
one-temperature accretion flow, for which / is significantly below unity, we can set q~ = where 
is given by Eq. (p^). The result can be expressed in terms of elliptical integrals, but to get a 
rough estimate of Qinc{R') we assume that R' <C R, i.e. we restrict our analysis to radii well inside 
the transition radius. Then to zeroth order in {R'/R) the ratio of the incident to the viscously 
dissipated energy evaluates to 

= 0.01 f) . (3-23) 

asTciR') \rtr) 

Thus, irradiation contributes at most a few percent to the energy output of the thin disk. This 
result was confirmed by integrating Eq. (3-22) numerically, with the correct expression for q~ 
in the case of a one-temperature flow. We conclude that irradiation is entirely negligible in this 
scenario, because the inner thin disk subtends a small solid angle from the point of view of the hot 
flow, and most of the energy is dissipated in the thin disk. By comparison, in the disk and corona 
model (e.g. Haardt &: Maraschi 1991), where irradiation is known to play an important role, the 
hot corona lies directly on top of the cold gas, so that roughly half of all photons emitted in the 
corona are absorbed by the disk. 

The expression for the energy incident on an outer thin disk differs from Eq. ( |3-22 ) only in 
the limits of the outermost integral; in this case we integrate R from the last stable orbit, 2>Rschw^ 
to Rtr- Again, we restrict ourselves to regions away from the transition radius, i.e. R' ^ R, and 
assume that q~ = q^ ■ Then to first order in R/R', the ratio of the incident to viscously dissipated 
energy per unit area is independent of R': 


Qinc{R') 

<tbTc{R') 

Numerical calculations show that even at R' ^ 


^ i?, the result increases by at most ~ 20%. 


(3-24) 


This calculation shows that irradiation of the outer thin disk becomes important only at large 
transition radii, where, according to our results from §3.3 the presence of the thin disk does not 
affect the hot flow. Note also that large transition radii are allowed only for low accretion rates. 


when the flow becomes strongly-advection dominated, and the result in Eq. (3-24) decreases by 
an additional factor of (1 — /). On the other hand, if the transition occurs within lO^Rschw, 
irradiation will increase the flux from the thin disk by at most a factor of 2. This has only a small 
effect on our calculations of the critical transition radius, and qualitatively our results remain the 


same. 


4. Discussion 
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In this paper we have explored how non-local radiation transfer affects the properties of hot 
one- and two-temperature advection-dominated solutions described by Esin et al. (1996), Narayan 
& Yi (1995) and Abramowicz et al. (1995). Since these solutions are optically thin, radiation 
emitted at one radius could in principle induce significant cooling or heating of the accreting 
gas at some other radius, through Compton scattering. We considered three different sources of 
non-local radiation; a luminous inner region of a hot disk {R < lOORschw)^ a thin disk inside a hot 
one-temperature gas, and a thin Shakura & Sunyaev disk surrounding a hot advection-dominated 
flow. 


We found that two-temperature solutions are not affected by the radiation from their hot 
interior. These solutions are strongly advection-dominated everywhere except near the outer 
boundary of the disk, Rcrit-, (see Fig. l[c]), and the inner part simply does not produce enough 
hot photons to induce significant heating of the outer layers. On the other hand, the gas density 
of the flow at ii > lOORschw is so small, that with electron temperatures of Tg ^ 10® K, Compton 
cooling by non-local radiation can not play an important role either. 

The inner parts of one-temperature flows are cooling-dominated and produce more radiation 
than their two-temperature counterparts. However, single-temperature solutions are limited to 
lower values of the mass accretion rate (see Figure l[a]), and since hot photons are produced 
mainly by Comptonization, which is not an efficient process at low M, no heating of the outer 
flows occurs. On the other hand, we find that Compton cooling by the photons from the interior 
dominates over local cooling processes for 10‘^Rschw ^ R ^ W^Rschw and M ^ MEdd- 
However, since the bulk of the observable radiation comes from inside lOORschw, the changes in 
the flow properties caused by extra cooling do not appreciably affect the overall spectrum of the 
system. 

Compton cooling by the radiation from the inner thin disk does have a significant impact on 
one-temperature flows, since it occurs at higher M, when Compton scattering is very effective. 
The extra cooling induced by the thin disk photons is so strong, that the maximum accretion rate 
for which the one-temperature model has an equilibrium solution decreases by a factor of ~ 10 
(Figure 3[a]). 

Finally, we found that the radiation from an outer thin disk does not have any effect on a hot 
flow that extends beyond ~ lO^Rschw At this distance, the gas is too cold for significant Compton 
cooling to take place. When the transition radius between the hot and cold gas moves closer to 
the accreting object, the non-local cooling of two-temperature gas near the boundary can become 
very large. However, this effect is limited to the layers of the hot flow very close to the transition 
between the hot and cold gas; away from this boundary, the ratio of non-local to local cooling 
decreases as Therefore, external cooling does not affect the regions of the hot gas where 

most of the observable radiation is produced and can be disregarded until it is so strong that it 
forces the hot flow near the boundary to cool to a thin disk configuration. When this happens, the 
transition radius decreases and the hot gas in the new boundary layer experiences even stronger 



- 19 - 


external cooling, which causes it to collapse in its turn. The resulting runaway process continues 
until the entire flow assumes the thin disk conflguration. This phenomenon can be important if 
the hot advection-dominated flows are formed by evaporation of the thin disk near the accreting 
black hole (Narayan & Yi 1995). In that case, the maximum accretion rate at which such hot flow 
can form is in fact less than Merit computed based on the local assumption. 


I would like to thank Ramesh Narayan and Eve Ostriker for many useful comments and 
discussions, and my referee for his/her friendly and helpful suggestions. This work was supported 
in part by NSF grant AST 9423209. 
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A. Modified Derivation of the Energy Eqnation And Adiabatic Exponent for the 

Accreting Gas 


One of the conditions that the physical parameters of the accretion flow must satisfy is the 
conservation of energy. Following Abramowicz et al. (1988) and Narayan &: Yi (1994) we can 
write this condition as 

where s is the entropy per unit mass of the gas and other quantities are defined in §2.1. 


(Al) 


To derive the expression for ds/dR we model the gas in the accreting flow as a combination 
of an ideal monatomic gas (a justified assumption since the accretion flow consists mainly of 
ionized hydrogen) and tangled magnetic fields. Radiation pressure was shown to be unimportant 
in one-temperature flows (Esin et al. 1996) and since two-temperature solutions are even less 
luminous, the radiation pressure can safely be ignored. Then the total pressure in the gas is the 
sum of the thermal ideal gas pressure, and the pressure due to the magnetic field (Eqs. |2-2| and 
12-31) and can be written as 


pkT 

p = Pg + Pm = -h — 

pm-u 247r 


pkT 

(5pmu 


(A2) 


where we have defined T = Ti, p = piPe/{Pe + PiTe/Ti) and (3 = Pg/p. In general, p is a function 
of both Ti and Tp, however, for simplicity we restrict ourselves to the two limiting cases when 
Tj > Te and p = pi, or R = Te and p = piPe/ipe + Pi). 

The internal energy of the gas is the sum of the kinetic energy of the particles (since the gas is 
assumed to be monatomic) and the energy stored in the magnetic field. The total internal energy 
per unit mass is then 


u = 


3 /cT 1R2 


3 kT 3(1-/3) kT (6-3/3) kT 


+ 


2 pniu p Svr 2 pruu (3 prUu 2/3 
In a quasistatic process, the first law of thermodynamics requires that 


pniu 


Tds = du + p dV = 


du 

&T 


dT + 


V 


du \ 

WJ: 


dV+pdV, 


(A3) 


(A4) 


where V = 1/p is the volume per unit mass. But tt is a function of T only, which means that the 
first term vanishes, {du/dV)T = 0. The second term we evaluate as 


du 

'Wf 


du 

It 


(6 - 3/3) k 


Einally, we divide both sides of Eq. 


Y ui 2/3 pmu 
1) by T and obtain the following relation for the entropy 


(A5) 


, (6 — 3/3) k dT 

ds = ^ ^ -^ 


pk 


dl- 


2(3 pniu T (3pmu \p 


(A6) 













- 21 - 


For an adiabatic process we set ds = 0. Following Clayton (1983, Eq. [2-121c]) we define the 
corresponding adiabatic exponent 7 as 


7 = Ts = 


I' 1 k 

V/3 fJ-rriu 




/ ( 6 - 3 / 3 ) k \ 
2/3 iimu} 


+ 1 = 


8-3/3 

6-3/3 


(A7) 


With this definition, we integrate Eq. 
to a constant factor: 


k 


to obtain an expression for the entropy of the gas np 
1 


s = 


\D.{(?gP^ '^) + const. 


(A8) 


Idpruu 7 - 1 

Now we are in a position to write down the final form for the energy conservation eqnation. 
Evalnating the derivative of s with respect to R and substitnting it into Eq. (^), yields the same 
expression as given by Narayan & Yi (1994) 


pv 


del 


/3(7 - 1) dR 

bnt with a different valne for 7, namely Eq. (A7). 


. 2 ^ 

dR 


-vcZ— = fq^ 


(A9) 
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B. Approximate Analytic Expression for the Viscous Heating and Local Cooling 


We need to obtain simplified expressions for the rate of viscous energy dissipation and local 
cooling of the accreting gas, as functions of the global parameters, m, rh, a, (3 and r only. 


B.l. Two-Temperature Accretion Disks 


In two-temperature flows the ion temperature is always nearly virial and scales as M/R oc 1/r. 
The electron and ion temperatures are equal when r ^ 1000, but Tg stays nearly constant at 
~ 10® K in the inner region of the disk. We find that the sum Tj -|- Tg is well approximated by 


Ti + T, 


10^^ K ( (3 

0^ 


(Bl) 


for all values of m, a, and /3. Substituting this expression into Eq. (^) gives us C3 ~ 0.3. 


Inverting Eq. (2-7), we obtain an expression for e' in terms of C3: 

9c3a^ 


, 1 f 9a^ X ' 

^ “ 2 1 “ “ 2 “ , 


X = 


(B2) 


The second term in brackets can be neglected as long as a ^ 0.5. Then Eq. (B2) simplifies to 

e' ~ 0 . 8 . 


Substituting these results into Eq. (2-5), we find that the total energy per unit volume 
dissipated in the disk can be written as 


~ 8.1 X 10^® —TT ^cm 


m 


-1_-3 


(B3) 


As required by the energy balance equation, the rate of local cooling is given by q~ = {3 — f)qt- 
In general, / is a function of r as well as m, a, and f3, but for simplicity we use the radially 
averaged value which can be easily estimated by integrating q~ over the volume of the disk and 
comparing the result with the total emitted luminosity, L 2 -t- 


/ OO 

qt 2 H 27 rRdR. 

- min 

We adopt the value for Ltot estimated by Mahadevan (1997) 

L2-t ^ 


(B4) 


1.3 X 1038 (^) (^) ergs-3, rh Z lO'^a^, 
2.5 X 10^^ mm ergs“^, m ^ lO^^a^. 


(B5) 
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Combining Eqs. (| 
the local cooling 


, and (B5) and setting = 3, we obtain the following expression for 


Q = 


(1 - f)qt = 9.3 X 10 

1.2 X 10^1^ 

m‘^Oi 

2.3 X 10^'^^ 


-18 ^ 2-1 


ergs ^cm 


m 


for rh ^ 10 
for rh ^ 10' -V. 


(B6) 


B.2. One-Temperature Accretion Disks 


Two-temperature flows become effectively single-temperature when r > 1000, since there 
Tg = Ti (Esin et al. 1996, Fig. l[b]). Moreover, in the region 30 ^ r ^ 1000 the gas temperature, 
T, in one-temperature disks is nearly as high as Tj in equivalent two-temperature solutions, and 
the values of the advection parameter / are similar as well. Therefore, in the outer regions Eqs. 


and (|B6D apply equally well to both models. 


Inside ~ 'i^Rschw, as one-temperature disks become thinner, denser and cooling-dominated, 
both and q~ increase faster with decreasing radius than the corresponding quantities in 
two-temperature disks. In this regime, we can use Eqs. (^) and (|B6|) as lower limits. 
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Figure Captions 

Figure 1. (a) Critical mass accretion rate, Merit, as a function of radius for one-temperature 
(dashed line) and two-temperature (solid line) models, plotted for two different values of the 
viscosity parameter a. Hot advection-dominated solutions exist below the respective lines; above 
the lines the only equilibrium configuration available to the accretion flow is a standard thin disk, 
(b) Radial temperature profile for a flow with a fixed accretion rate of MEdd- (c) Contours of 
constant advection parameter / for a two-temperature accretion flow. The critical accretion rate 
is shown as a heavy line, (d) Contours of constant advection parameter / for a one-temperature 
accretion flow. The critical accretion rate is shown as a heavy line. 


Figure 2. Critical mass accretion rate (heavy line). Merit, plotted in Eddington units, as a 
function of radius for two-temperature (a, c) and one-temperature (b, d) models, (a) and (b) 
Above the dashed line lies the region where significant Compton heating by the photons emitted 
at i? < WORschw is expected to occur, based on an analytic approximation (Eqs. |3-5| and 
ll3-7[ ). The dotted lines are the contours of constant qni/q~ computed numerically, labeled by the 
values of this ratio. Positive numbers correspond to a net heating, negative numbers indicate 
that non-local radiation produces a net cooling. Note that the ratio of non-local to local heating 
is always significantly below unity, (c) and (d) Shaded area corresponds to the region of the 
accretion flow where local cooling processes are less efficient than Compton cooling by the photons 
emitted at i? < WORsehw, computed numerically. The dashed lines are the contours of q~i/q~ = 1, 
estimated analytically (Eqs. ||3-9| and |1-10[ ). 


Eigure 3. (a) Critical mass accretion rate, Merit, plotted in Eddington units (solid heavy line), 
for a one-temperature model with a IOMq black hole, viscosity parameter a = 0.1, and ratio of 
gas pressure to total pressure (3 = 0.5. For accretion rate values above Merit the only equilibrium 
configuration available to the accretion flow is a standard thin disk. The region of the hot accretion 
flow for which non-local cooling due to scattering of the cold photons emitted by the thin disk is 
more efficient than local cooling processes is shaded. When the non-local cooling is included in 
the calculations, the critical accretion rate decreases. The new Merit curve is plotted as a dashed 
heavy line, (b) Compares the temperature profile of the flow with (thin line) and without (heavy 
line) extra cooling due to the inner thin disk. 


Figure 4. The ratio of non-local cooling rate due to scattering of the radiation emitted in the outer 
thin disk, to the local cooling rate in the gas near the transition layer between the hot and cold 
flows, plotted as a function of the transition radius for two-temperature (a) and one-temperature 
(b) models. Different curves are labeled by the values of the accretion rate in Eddington units. 
Note that for R^r ^ lOORsehw, this ratio is significantly smaller in one-temperature models which 
have much stronger local cooling than their two-temperature counterparts. 
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Figure 5. Critical mass accretion rate, Merit, plotted in Eddington units (solid lines) for 
two-temperature models with viscosity parameter a = 0.1 and 0.3. The dashed lines correspond to 
the minimum allowed values of the transition radius between the hot flow and the outer thin disk. 
If the transition occurs in the region to the left of the dashed line, Compton cooling of the hot gas 
by the thin disk photons is so strong that the entire flow settles to the thin disk configuration. 
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